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SUMMARY 

A preliminary theoretical analysis was conducted of expansion- tube and expansion- 
tunnel simulation whereby the postnormal shock stagnation conditions of a model were 
matched to those of flight for both a perfect gas and equilibrium real air. The selected 
simulation mode permitted large variations from flight conditions for the test- section 
ambient temperature and pressure but only a small excursion for ambient velocity or 
density. 

A brief treatment of the perfect- gas expansion tunnel with simulation demonstrated 
both gains and penalties in driver performance, depending on the Mach number, nozzle 
area ratio, driver gas, and simulation temperature ratio. 

A more exhaustive analysis of an equilibrium real-air expansion tube employing 
simulation with ambient-gas temperatures of 1000° and 2000° K was completed for equiv- 
alent flight velocities from 20 000 to 50 000 feet per second (6.096 to 15.24 km/sec) and 
altitudes from 100 000 to 300 000 feet (30.48 to 91.44 km). The primary advantages of 
simulation were found to be large increases in testing time, initial test-gas slug length, 
and velocity -altitude performance. The primary disadvantages were a small degree of 
test- section dissociation and an increased peak cycle dissociation for the highest flow 
velocities and altitudes. 

For equilibrium flow the gas state behind a shock wave was found to be closely 
simulated by all strong shock waves and to be well simulated even by weak waves with 
flow deflections as low as 30°. 

INTRODUCTION 

The characteristics of expansion tubes and expansion tunnels operating to generate 
test-gas flows wherein the test medium exactly duplicates flight conditions have been 
reported in references 1 to 4. Certain advantages and disadvantages are obvious, and the 
question naturally arises as to means for reducing particular disadvantages without either 
a concurrent unacceptable reduction in advantages or the promotion of additional 





disadvantages. The short testing time and the short initial extent of the test gas slug 
adjacent to the secondary diaphragm are two prime disadvantages of the expansion tube. 
One means of alleviating these particular drawbacks is the expansion tunnel wherein a 
nozzle is added to the downstream end of an expansion tube (refs. 3 and 4). 

Another attack on the same problems is the use of simulation in which the test gas 
does not exactly duplicate all the flight conditions; instead, the test-gas parameters are 
so chosen as to duplicate the most important features of the flight system under study. 

A typical application of simulation to a different problem is the use of shock tubes to 
study the blunt-body stagnation-point heating rate of ballistic missiles (ref. 5). By cor- 
rect duplication of the gas state at the stagnation point, the flow in a shock tube at a 
Mach number of 3 was employed to study experimentally a real-gas atmospheric reentry 
at hypersonic Mach numbers. 

A study of the application of simulation to expansion tube operation was begun at 
Langley Research Center in 1961. Since the more recent atmospheric tables (ref. 6 ) and 
thermodynamic properties of air (ref. 7) were not available in 1961, the 1959 ARDC 
atmosphere (ref. 8 ) and air properties of references 9 and 10 form the basis for the 
analysis. Some of the results of this analysis together with a brief perfect-gas study of 
simulation in an expansion tunnel are presented in this report. A real-gas analysis of 
the expansion tunnel both with and without simulation is also under way at Langley but 
results are not sufficiently complete to be reported herein. 

SYMBOLS 

A area ratio of exit to inlet of nozzle in expansion tunnel 

A 4 cross-section area of driver 

a speed of sound 

c v specific heat at constant volume 

E driver energy 

H total enthalpy 

h local enthalpy 

l d>^ Sj>^ S 2 lengths of driver, intermediate (driven) chamber, and accelerating 

(expansion) chamber, respectively 


2 


Mach number 


local pressure 


standard atmospheric pressure 


body nose radius 


ratio of entropy S to gas constant R 


local conditions at stagnation point of blunt body 


local temperature 


ideal test time (fig. 19) 


primary shock velocity 


fluid or flight velocity 


reciprocal molecular weight ratio, 


Undissociated molecular weight 


Dissociated molecular weight 


initial length of test-gas slug in intermediate chamber 


_ ri(M e ) rr 2 (M f )' 


(M f )[T 2 (M e ) 


r 2 (Mi) = i + Mi 2 


ratio of specific heats 


normal- shock detachment distance 


flow- deflection angle 


3 



reciprocal of normal- shock density ratio, Density before shock wave 

Density after shock wave 

ideal driver efficiency (eq. (28)) 
nondimensional time parameters, see reference 11 
viscosity 
local density 
shock-wave angle 

Subscripts: 

b local conditions behind normal shock 

e entrance to expansion-tunnel nozzle 

f exit from expansion-tunnel nozzle 

i denotes gas in state i in cycle, identified by circled numbers in figure 1 

°° denotes gas in atmospheric flight (free stream) 

( )i subscript signifies quantity to be evaluated in ith state 


e 

^o 

M 

P 

a 


SIMULATION THEORY 

Criteria for simulation are initially derived by matching the adiabatic -flow 
stagnation- point parameters behind a normal shock for a simulated free stream identified 
in figure 1 as region (5) to those of flight which are identified by the subscript (See 
following sketch.) 
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( 1 ) 


Since H5 = at the stagnation point, conservation of energy requires that 




h 5 


which may be manipulated to determine the simulation velocity ratio 



Conservation of momentum across the shock wave is expressed as 


Poc + P 


=o U co =P 


, + p, ur 

Doc Doo Doo 


( 2 ) 

(3) 


Combining this equation and the continuity equation with the following incompress- 
ible Bernoulli equation 


<v . + 1 vL * Pb. 


(4) 


which is assumed valid between the shock and stagnation point, yields the approximations 
for the stagnation pressure, as follows: 


Poo + P 


oo U oo 



(5) 


Ps 5 ~P 5 + P 5 u 5 2 ( 1 " x) (5a) 

Consequently, for equal stagnation pressures for both flight and simulation ^p Soo = PS5) 
the dynamic-pressure ratio is found from equation (5) to be 


The simulation density ratio 


£5 

Poo 


P5 U 5 
P M 


»-T + 


PoqUoo 


6 5 

T 


P5 


P 5 U 5 


P5/P00 re suits from combining equations (2) and (6) 
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( 6 ) 
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If the assumption of equilibrium and irreversible flow 
foregoing equations have a unique solution for any state °°, 
ships exist: 

For equilibrium composition, 

p 5 = p 5K h 5) W 

For equal total enthalpies (eq. (2)), 

u 5 = u 5( h 5) ( 9) 

For equilibrium normal shock, 

6 5 = 6 5(P5> h 5> u 5) ( 10 ) 

= e 5(P5’ h 5) ( 10a ) 

According to equation (7), 

P 5 = P5( e 5’ p 5’ h 5) ( n ) 

Thus, from equations (8), (10a), and (11), 

P5 = P5( h 5) ( lla ) 

For a dissociating and ionizing gas, however, the solution is iterative rather than 
in closed form. Consequently, the sections to follow discuss the easily obtained but 
approximate solutions which are sufficiently accurate for the high velocities (Mach num- 
bers) of interest. 

PERFECT-GAS SIMULATION FOR EXPANSION TUBE AND EXPANSION TUNNEL 


Sb 5 „ S 5 


R 


s — J is applied, the 


since the following relation - 


Since the purpose in illustrating the perfect-gas simulation is to indicate general 
performance trends at moderate values of the hypersonic restriction of 

Moo » 1 and M 5 » 1 is introduced. Thus, equations (2) and (7) reduce to 


U5 

Uoo 


1 


( 12 ) 



The other simulation conditions are then 


(13) 


p 5 _ T 5 

Poo Too 


( 14 ) 
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( 15 ) 


m 5 _ aoo _./T^ 

Moo a 5 jfTg 

' Consequently the prediction of the performance of an expansion tube or tunnel sim- 
ulating the conditions M^p^ is simply reduced to the computation of a tube or tunnel 

with Me = and Pe = p (■ — \ In particular, consider the ratio of the driver 

JTg 0 ^Too/ 

pressure to static pressure P4/P00- The governing equation is equation (36) of refer- 
ence 3 which is rewritten as 



Figure 1 of the present report gives the notation used for the expansion- tube cycle, and 
figure 1 of reference 3 gives that for an expansion tunnel. (Note that test state f of 
ref. 3 is replaced herein by state 5.) The driver-pressure ratio is then 



In figure 2 curves are drawn of the ratio (p^/p^) to (P4/P°o) 

T^Too T5=T 0O 

equation (17) with y 4 = 5/3, y = 1.4, ^ when = 10 or 25 


by using 
and 


Tg/Too = 1/3, 10/3, or 20/3. This figure is the perfect-gas analogy to a heated helium- 
driven facility using driver temperatures T4 of 3000° K and 7500° K with test- section 
air temperatures Tg of 100° K, 1000° K, and 2000° K to simulate flight in an atmos- 
phere where = 300° K. 

An approximate equation relating P^j/P^ to M^, Tg/T,*,, and A for Mg » 1 
and M*, » 1 is obtained by using equation (17b) of reference 3 for j3ri(Mgj: 
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over 


The variation of the ratio of the driver pressure (P4) T ^ T to (P4) T _ T 
a wide range on both sides of unity may be understood by consideration of the terms in 

equation (18). For small values of A the term 

of Moo and the driver-pressure requirements are reduced by a simulation in which 
T5 > T^. As Moo increases the relative contribution of the temperature- ratio term is 
reduced, and the driver-pressure ratio decreases further. However, as M^ increases 
still further, the bracketed denominator becomes smaller so that the driver-pressure 
ratio begins to rise. (Note that there is a limiting M^ when the pressure ratio P4/PC0 
becomes infinite.) 


7 tT~J dominates in the lower range 


As A increases with other conditions fixed, the driver-pressure ratio increases 
because the numerator of the fraction (eq. (18)) increases while the denominator 
decreases; and this pressure increase is further undesirably augmented by simulation 
when T5 > T^. Thus, in regard to the driver-pressure requirements there is a 
restricted range where the use of simulation with T5 > is advantageous. When 
T5 < Tco the ratio of (P4),p ^ to (p^^ _,p is greater than unity for 10 < Moo < 50 

and 1 S A ^ 10^ (fig. 2(a)). Consequently, on a driver-pressure basis, simulation with 
Tg > Too appears more advantageous in this Mach number and area range. 

Equilibrium- air calculations for simulation with A > 1 have beetf under way for 
some time at Langley Research Center, but are not considered in this analysis. The 
next section treats the particular case of Tg > T w and A = 1 for the expansiorijube. 


EQUILIBRIUM AIR SIMULATION FOR EXPANSION TUBE 


Conditions in Various States of Cycle 

Conditions in test-gas region 5. - To simplify the determination of ug/u,*, and 
P5/P00 f° r equilibrium real air, the following approximate scheme is used. First, the 
right-hand side of equation (6) is assumed to be unity. Such an assumption introduces 

only a small error for large u M since €<*, ~ e 5 « 1 and ~ — 9 « 1. Second, 

pu^ yM^ 

in equations (1), (2), or (7) the term 2h 00 y(i 00 2 W as assumed to be negligible compared 

to unity /i.e., M,*^ » — ; but the term 2hg/u 00 2 was re tained so as to increase the 

' ^ ’ 9 T 2 

accuracy of the solutions since Mg* is of the order of — & . Because the inde- 

5 

pendent parameters selected were p^, h^, Uoo, and Tg (rather than hgj, the reten- 
tion of the hg term required an iterative solution which, however, converged very 
rapidly so that acceptable accuracy was obtained after a single iteration to the following 
modified equations: 
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(19) 

( 20 ) 


h 5 ~ h 5( p 5’ T 5) ( 21 ) 

p 5 = p 5( p 5’ T 5) ( 22 ) 

Equations (13) to (22) were solved for the simulation state 5 when T5 = 1000° K 
and 2000° K. The ambient conditions simulated were for velocities u w between 20 000 
and 50 000 feet per second (6.096 and 15.24 km/sec) at densities and temperatures 
corresponding to flight altitudes of 100 000, 150 000, 200 000, 250 000, and 
300 000 feet. Since this calculation was performed prior to the release of the U.S. 1962 
standard atmosphere (ref. 6), values of the ARDC 1959 model atmosphere (ref. 8) were 
employed for flight conditions «>. The differences between the two atmospheres are 
significant only for the highest of the altitudes considered herein. 

The normalized simulation conditions of the test region (state 5) appear in fig- 
ures 3 to 7. The normalized velocity ratio (fig. 3) has an increased departure 

from unity for lower velocities and increased T5; this effect is a direct result of the 
retention of the hg/u^ 2 term which increases as Tg increases or u^ 2 decreases. 

For Tc = 1000° K, 0.975 < ^ < 0.997 and for T 5 = 2000° K, 0.922 < ^- < 0.991 in 

the velocity range 20 000 to 50 000 feet per second (6.1 to 15.24 km/sec). The normalized 
density ratio (fig. 4) has twice the variation of the velocity ratio (eq. (20)) with a maximum 
value of 1.05 and 1.18 for Tg = 1000° K and 2000° K, respectively, at the lowest velocity. 
The reciprocal molecular weight ratio Z5 (fig. 5), which (because the Van der Waal 
forces are negligible in the test region 5) is equivalent to the compressibility factor, is 
essentially unity for Tg = 1000° K. For Tg = 2000° K, Zg varies with altitude from 
1.00 up to 1.04. Note that figure 5 does not indicate a variation of Zg with u B at a 
given altitude. Actually, there is a small variation with ug, but except at the highest 
altitude, it is not discernible within the resolution of this figure, and hence an average 
value is plotted. 

The fraction of the total enthalpy due to ambient dissociation is less than 1 percent 
for altitudes below 250 000 feet (76.2 km) and Tg = 2000° K. However, at 300 000 feet 
(91.44 km) the dissociation energy for Tg = 2000° K is 3 percent and 1 percent of the 
total enthalpy for Uoo = 20 000 and 50 000 ft/sec (6.096 and 15.24 km/sec), respectively. 
Such small percentages should not distort the resultant model flow fields significantly. 
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The normalized pressure ratio P^/Pc-, (fig- 6) exhibits a larger variation with 
altitude for given velocity than either the density or velocity. This variation arises 
mainly because T^ is a function of altitude, but it is also influenced slightly by Z5 

(for Tg= 2000°K). Consequently, since — = Z 5 the variation of Pg/p^ with 

altitude will be larger than that of Pg/Poo* 

The normalized viscosity ratio M5/M00 (fig- V is also dependent on altitude. 

Since the Reynolds number ratio ^ -- U - ! - ^°° = — the viscosity ratio essentially 

p 5 Poo U °o u 5 M5 

determines the Reynolds number ratio. Figure 7 thus indicates a lowering of the 
Reynolds number by factors up to 6, which occur in simulation at 2000° K. Since wave 
attenuation, boundary- layer growth, and test-gas velocity profile are all influenced by 
Reynolds number, this lower Reynolds number is an undesirable feature of simulation. 

(See ref. 3.) 

The dimensionless entropy Sg/R is plotted in figure 8 as a function of altitude. 
Although only a single line is shown for both T5 = T^ and Tg = 1000° K, the calcula- 
tions for T5 = 1000° K actually plot a family of curves, which are not separable within 
the definition of the figure. The two extremes of this family are shown for Tg = 2000° K; 
the difference in Sg/R between u^ = 20 000 ft/sec (6.096 km/sec) and 
Uoo = 50 000 ft/sec (15.24 km/sec) is very small. 

Conditions for region 2 .- The conditions for region 2 were determined from those 
in region 5 by using the charts of reference 11 for equilibrium one- dimensional unsteady 
isentropic flows, together with the values of ug and Sg/R from figures 3 and 8. In 
certain cases a slight extrapolation of the data in references 9 and 10 was necessary to 
obtain the conditions in region 2 but the errors so introduced are believed small. No 
computations were completed for u > 35 000 ft/sec (10.668 km/sec) at an altitude of 
100 000 feet (30.48 km) because the required extrapolations became too large. Since 
these computations were completed, reference 7 has been published and the data therein 
would eliminate the necessity for such extrapolation. 

Figure 9 shows that the flow velocity U2 increases as Tg increases for a given 
u w . This trend is to be expected since the total enthalpy multiplication ratio H5/H2 is 
reduced by expanding only to Mg instead of Moo- (See eq. (16) and fig. 7 of refer- 
ence 1.) Note, however, that Mg of this reference is equivalent to the Mg herein 
only for Tg = T M . Consequently, since H2 ~ U£^, a higher value of U2 is required to 
produce the same Hg. 

Temperatures T2 for simulation are also generally increased for T5 = 1000° K 
and 2000° K (fig. 10). The explanation of this variation follows along the same lines as 
that in the previous paragraph since H2 ~ 21^, and in general for small dissociation T2 
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varies as h^. However, for appreciable dissociation the temperature is depressed. 

This depression is responsible for the approximate invariance in T2 with Tg at the 
higher velocities. 

In contrast to both the temperature and velocities in region 2 and pressures in 
region 5, the pressures for simulation in region 2 with T5 > T,^ are lower for a given 
velocity and altitude. These pressures, normalized by standard atmospheric pressure 
Pq in figure 11, are approximately a decade lower when T5 = 2000° K than when 
T5 = T^. Thus, the increased simulation entropy more than compensates for the 
increased T2 with the net result that pg decreases. 

One of the principal disadvantages of simulation with Tg > T^ is evident from 
the data shown in figure 12; namely, the appreciable increase in Z2 for Tg = 1000° K 
or 2000° K. When Tg = T M , Z2 = 1.00 at u,*, = 20 000 feet per second 
(6.096 km/sec), whereas when T5 = 2000° K the values of Z2 vary from 1.06 to 1.15. 
Similarly, at u*, = 50 000 feet per second (15.24 km/sec), 1.09 = Z2 = 1.18 for 
T5 = Too, whereas 1.20 S Z2 = 1.33 for T5 = 2000° K. Since by the definition of Z 
used herein, the mass fraction dissociated is essentially Z - 1, the dissociation in 
region 2 has been approximately doubled in a change from Tg = T^ to Tg = 2000° K 
at u^ = 50 000 feet per second (15.24 km/sec). If these atoms are not recombined 
(i.e., if flow "freezes") during or following the unsteady expansion, then a larger degree 
of dissociation will exist at the test section than is indicated in figure 4. The fact that 
Z2 is higher and p£ is lower for simulation than for duplication increases the possi- 
bility of a greater degree of nonequilibrium or even of a frozen test gas for simulated 
testing as compared with that for duplication. Countering these adverse factors and 
tending to maintain equilibrium for simulation is the much slower rate of expansion which 
the test gas experiences because the test time is longer and the last part of the test-gas 
slug therefore enters the unsteady centered expansion wave farther from the singularity 
at the origin. The test time and slug length are discussed in a subsequent section of this 
paper. Of course, even the values of Z2 a t 2000° K are much less than the maximum 
Z occurring before the nozzle expansion in a reflected shock tunnel. For example, at 
u^, = 30 000 feet per second (9.144 km/sec) and an altitude of 200 000 feet (60.96 km) a 
reflected shock tunnel would require a Z at the nozzle entrance of approximately 1.35 
as contrasted to Z 2 * 1.16 for Tg = 2000° K. Thus for u,*, = 30 000 feet per second 
(9.144 km/sec) at 200 000 feet (60.96 km) altitude the shock-tunnel dissociation would 
still be about twice the expansion- tube value with 2000° K simulation. 

Conditions for region 1. - The shock velocity Ug^ and pressure pj were deter- 
mined from the charts of reference 11 for Tj = 300° K. The trend of Us2 for simula- 
tion (fig. 13) is, as expected, similar to that of U2 (fig. 9) with Usj increasing as Tg 
increases. Also, p./pQ (fig- 14) decreases in approximately the same manner as 
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P2/ p 0 

vacuum 


(fig. 11). The fact that lower values of Pj/Pq are required means an increased 
capability for the intermediate chamber with minimum pressures of Pi/Pq of 


the order of 10"® for T5 = 2000° K as compared with values on the order of 10"^ for 
T 5 = T OO (fig. 14). However, such values of Pj/Pq are fairly easily obtained and are 
significantly above the values of Pjq/Pq (which are the minimum pressures in the over- 
all apparatus) to be discussed in the following paragraph. 


Conditions for region 10 .- Perfect helium was arbitrarily assumed for both 
regions 20 and 10, and the approximation of a strong shock was invoked. Consequently, 
equation (3) of reference 1 was used to obtain P2o/ p 10’ an< * p lo/ p O then f°^ ows since 

p 20 = p 5 = (j^) p °o' In ^ his case ’ va * ues °f t 5 > Too are beneficial because Pjo/ p O 


increases as T5 increases (fig. 15) and reduces the vacuum requirements of the expan- 
sion chamber. This increased pressure is a substantial advantage since the difficulty in 
obtaining Pio/Po = 3 X 10 _ ® is much less than in obtaining P^o/PO = 2 x 10“^. (These 
values are for T5 = 2000° K and T5 = Too, respectively, at Uoo = 50 000 feet per 
second (15.24 km/sec) for an altitude of 300 000 feet (91.44 km). The assumptions of 
perfect helium and a strong shock are naturally questionable at these low pressures. 


Conditions for driver region 4 .- Three driver gases are considered to illustrate 
conditions in the driver region 4: 

(1) A combustion driver consisting of a mixture of 75 percent helium with 25 per- 
cent of a stoichiometric mixture of oxygen and hydrogen (fig. 16); 

(2) Hydrogen heated to T4 = 833° K (fig. 17); 

(3) Helium heated to T4 = 4000° K (fig. 18) 

These drivers might be considered typical of those heated, respectively, by combustion 
(ref. 12), electric resistance heaters (either external or internal, refs. 13 and 14), and 
electric -arc heating to moderate T4. Appreciable reductions in driver-pressure ratio 
p 4/ p 0 resu lt from simulation with a decrease of approximately a decade apparent for 
the three driver gases. The driver pressures required for given T5, u,*,, and altitude 
are nearly equal for the combustion and hydrogen drivers (figs. 16 and 17), whereas pres- 
sures for the heated helium are somewhat lower (fig. 18). 


Lengths of Component Sections 

A wave diagram of the component sections in the expansion tube is illustrated in 
figure 19. The ratio of the intermediate chamber length l to the accelerating cham- 
ber length l §2 is chosen so that the reflected wave originating from the intersection 
of the upstream edge of the expansion fan with the entropy discontinuity (point III of 
fig. 19) arrives at the test section concurrently with the downstream edge of the expansion 
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fan (point n). From this wave diagram the time interval between the secondary dia- 
phragm rupture and point HI may be simply related to the length Zg-p the conditions in 
region 2, and the primary shock speed. Similarly, the time interval to point II is 


l 


S 2 


„ „ . Since these time intervals 

u 5 " a 5 

ref. 11 ), it follows that the ratio of 


may also be expressed as f ( 77 ) (see eq. (15) of 
?Sj to ^2 may be expressed as 


l Jh 

l S2 


Uc 


u 5 “ a S 


a 2 

r a 5 

u Si " u 2 

L a 2 


-. 1/2 


(23) 


Charts of rj from reference 11 were used to evaluate equation (23) and the resulting 
ratios l g^yz gg are shown in figure 20. Increased intermediate chamber length ratio 
l Sl /l S2 is required for Tg > Too with nearly a decade increase required for 

Tg = 2000° K over values for Tg = T^. However, this increase is not of major con- 
cern for practical designs except at the lower velocities; for example, for Tg = 2000° K, 
Z g ^jl g 2 ~ 0.3 at Uoo = 20 000 feet per second (6.096 km/sec) and Z g ^jl g 2 = 0.02 
at Uoo = 50 000 feet per second (15.24 km/ sec). 


The values of Z Sj/^ S 2 °* 20 for Tg = T^ are not identical with those of 

reference 15 (fig. 11) or with the ratios obtained from figures 14 and 15 of reference 1 
because of the approximate integration scheme used in those references. However, the 
variation is small except for one case where an error in integration resulted in an incor- 
rect value in references 1 and 15 for Z Sl at an altitude of 250 000 feet (76.2 km). 


The driver chamber length Z q was chosen to match the intermediate chamber 
length Z gj so that the initial reflection from the end of the driver would pass through 
point HI, figure 19. The governing equation, which follows, was derived by equating the 


time necessary for the entropy discontinuity to travel the distance 



to the sum 


of the times required for the reflected wave to: (a) traverse the driver in a negative 
direction, (b) traverse the driver expansion fan, and (c) traverse region 3 to point HI. 

y 4 +i 

iD L % - “2\ a 4 (J^> 

*S X \ a 2 /2U Sl \ a 4/ 


(24) 


The sound speed ratio for the driver expansion is 


_ j _ y 4 ~ 1 ^2 
a 4 2 a 4 


(25) 
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The length ratios Ijyjl g^ are shown in figures 21 to 23 for the three driver gases of 

figures 16 to 18. In all cases the ratio l d/^Sj decreases as T 5 increases. At 
T 5 = 2000° K the driver length ratio l d/^Sj approximately one-half that requirecf 
for T 5 = T OO* However, this decrease in Ijyjl g^ is not sufficient to offset the increase 
in isi/^S 2 > w ^h the ne t res ult that the ratio is an order of magnitude larger 

for T 5 = 2000° K than for T 5 = T^. The significance of the ratio *s dis- 

cussed further in a subsequent section. 

Test Time Per Unit Length 

One of the most critical factors in expansion- tube operation is the length of the 
ideal test time, which is indicated as At 2 (fig. 19). The ideal test time per unit accel- 
erating chamber length S 2 plotted in figure 24. Significant gains are realizable 

for increasing T 5 ; at T 5 = 2000° K, Atjl^ I s roughly larger by a factor of 3 than 
that for T 5 = Too at the same u w and altitude. In a practical application, gains larger 
than those just mentioned may be possible since the actual test time will be less than the 
ideal time by the time required to establish the flow about a model once the fluid in 
state 5 reaches the model, as well as by the time lost due to mixing and other effects. 
These flow- establishment and mixing-loss times should not change appreciably with 
simulation. Thus, consider a case where the actual test time was only 25 percent of the 
ideal test time when T 5 = T*,. Then, for the same lost time and an ideal test time with 
simulation larger by a factor of 3, the actual simulation test time would be 900 percent of 
that for T 5 = T^. 

Although simulation increases At 2 ^ $ 2 ’ the opposite trend results for At 2/7 Sj 
(fig. 25). Consequently for given flight conditions °° and ideal test time, a longer inter- 
mediate chamber is required. The values of At/Zg^ for T 5 = Too are roughly 5 to 
7 times larger than those for T 5 = 2000° K. 


Initial Test Slug Length 


The criticalness of the secondary diaphragm rupture because of the very limited 
axial distance occupied by the test-gas slug while in regions 1 and 2 has been pointed out 
in references 1, 3, and 4. The quantity z is defined as the initial extent of the test gas 
slug in region 1 ; that is, 


P 5 Ac 

z = — — -T— Uc At 

Pi A 1 5 


( 26 ) 


where for purposes of this analysis At = At 2 in equation (26). The normalized ratio 
z / 1 Sg exhibits very large gains for simulation with T 5 > T^ (fig. 26). In fact, for 
T 5 = 2000° K, z jl g^ is approximately two orders of magnitude greater than that for 
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T 5 = Too. Of course, even with such gains, the absolute value of z/Zgg is very 
small at the higher velocities. The possibility of further increasing z jl g 2 by a com- 
bination of simulation with expansion tunnel operation is a fruitful field for investigation 
since both modes of operation individually increase z/Zg 2 significantly. 


Driver Considerations 


Arc driver energy .- Although for the combustion or resistance heated driver gases 
the states in region 4 available to the experimenter are primarily determined by p^, 
which is in turn limited by the strength of the driver vessel, the user of an electric- arc 
discharge-heated driver gas is additionally limited by the amount of energy E in his 
condenser bank or inductive field (assuming 100 percent arc efficiency). A convenient 
parameter to evaluate this energy restriction is the driver energy density E /A 4 Z p, 
which for a perfect driver gas is 


E 

a 4^D 



(27) 


Since the driver energy density for a perfect gas is simply related to the driver 
pressure p^, the variation in energy density of figure 27 for a helium driver at 
T 4 = 4000° K follows the pressure variation of figure 18. Again for otherwise compar- 
able conditions a decrease in E /A4.I D of approximately a decade occurs between 
T 5 = T^ and T 5 = 2000° K. This energy parameter will vary with T 4 . (See, for 
example, the variation of p 4 with T 4 in fig. 12(a) of ref. 1.) Even with simulation, 
large values of E are required for expansion tubes of practical size since for 

T 5 = 2000O k, Uoo = 50 000 ft/sec (15.24 km), * 3 x 10 4 joules/ft 3 

■**4^ D 

(l X 10® joules/m 3 ) at 300 000 feet (91.44 km) altitude and 1 x 10® joules/ft 3 
(4 x 10® joules /m 3 ) at 100 000 feet (30.48 km). 

Although the energy density parameter permits the matching of the energy source 
E with the driver dimensions, this parameter gives no indication of the gains or losses 
from simulation, since the driver length l D required also varies with the degree of 
simulation. However, such an indicative parameter is the ideal driver efficiency r? 0 , 
which is the ratio of the kinetic energy of the test- section gas to the driver energy E: 

_ (K)(P5»5At2A 5 ) 

'O E 


From figure 28 it is evident that simulation does improve the ideal efficiency with 
typical increases of half a decade at the higher velocities. Caution should be exercised 
in the use of this figure because of the inherent assumption that the lengths l j) and 
1 can be continuously varied to obtain the previously discussed optimum ratios of 
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I j -yjl and l g ^Jl gg. However, even with the moderate flexibility obtainable by the use 

of a number of finite sections to comprise the various chambers, most of the gains due to 
an infinite flexibility should be attainable since configurations both with and without 
simulation would be penalized. 

Another factor to be weighed is the ideal length of the driver required as compared 
with the practical physical size from consideration of viscous effects, diaphragm rupture, 
and other factors. For example, the ratio l -q/1 S 2 formed from figures 20 and 23 is 
plotted in figure 29. Note that larger values of Ijy/l §2 are required for simulation; 

for example, at 50 000 feet per second (15.24 km/sec) —2 * 1.5 X 10"^ at T 5 = T,*, 

l ?S2 

and -5- * 2 X 1(T® at T 5 = 2000° K. Assuming a value of l go = 100 feet (30 m) and 
Z S 2 

a diameter of 1 foot (0.3 m), then the ideal ~ 0.015 foot (0.005 m) for T 5 = Too and 
0.2 foot (0.06 m) for T 5 * 2000° K. Such lengths appear unrealistic for a 1-foot- 
diameter driver; in fact, a driver length- to- diameter ratio of unity might be a practical 
minimum. In such a case the ideal efficiencies of figure 27 would be reduced by factors 
of 0.015 and 0.2, respectively, so that the half-decade ideal efficiency gains for simula- 
tion become "practical efficiency" gains of a decade and a half. 

Typical maximum-velocity-performance comparison .- The maximum Uoo per- 
formance of an expansion tube with and without simulation was determined from fig- 
ures 16 and 17 for combustion (75 percent He) and hydrogen (T 4 = 833° k) driver with a 
maximum P4 = 20 000 psi (l.36 X 10 2 atm) and also from figure 27 for a helium driver 
(T 4 = 4000° K) with E = 2.5 x 10® joules, A 4 = 1 ft^ (o.09 m 2 ), = 1 ft (0.3 m). 

Note that these three drivers give almost identical performance for T 5 = T^ (fig. 30). 
However, as T 5 increases the gains in u,*, are smaller for the combustion and hydro- 
gen drivers than for the helium driver. For T 5 = 2000° K a general rule of thumb 
describing the simulation gains for all three drivers might be an increase in of 

roughly 10 000 feet per second (3.048 km/sec) at the same altitude, or reduction in alti- 
tude of 75 000 feet (22.86 km) at the same u^. Such performance increases are very 
significant gains in facility capability. 

EFFECT OF SIMULATION ON FLOW FIELDS 
Equilibrium Case 

The simulation method investigated was selected to duplicate equilibrium condi- 
tions at the stagnation point of a blunt body in hypersonic flow. Since the gas states do 
not change appreciably in the stagnation region bounded by the normal shock and the 
blunt-body stagnation point, the conditions behind the normal shock should also be closely 
duplicated. For a given body nose radius r, the shock detachment distance 6 is 
inversely proportional to the stagnation- region density ratio; consequently, the ratio 
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(6/r) 5 

j-j ■— will be proportional to the ratio P^jP <*>• The data of figure 4 imply that the 

ratio 65/8^, at T & = 1000° K decreases from 1.05 at 'J M = 20 000 ft/sec 

(6.096 km/sec) to 1.01 at u^ = 50 000 ft/sec (15.240 km/ sec). At T 5 * 2000° K the 

discrepancy is increased, with a variation in 65/6^ from approximately 1.15 at 

Uoo = 20 000 ft/sec to 1.02 at u w = 50 000 ft/sec. If the maximum deviation acceptable 

was 10 percent, then simulation at Tg = 2000° K would be valid only for 

Uoo > 25 000 ft/sec (7.620 km/sec). 


Simulation decreases the test section Mach number M5 approximately as the 
square root of T00/T5 (fig. 31). Such a drastic reduction in Mach number might be 
important if the reduced Mach number was not still hypersonic and if the test fluid did 
not satisfy the simulation requirements herein. However, even for the lowest velocity 
considered in this paper, the Mach number is above 7 (u^ = 20 000 ft/sec (6.096 km/sec), 
T5 = 2000° K). References 16 and 17 indicate that for strong shock waves in a real gas 
the total enthalpy is the most important parameter influencing shock shape and total 
enthalpy is conserved in simulation. 


Shock-wave polars for the simulated flows were computed by the equilibrium 
normal-shock program of reference 18. The results are shown in figures 32 to 35 where 
the left- and right-hand ordinates refer to Uoo = 20 000 and 50 000 ft/sec (6.096 and 
15.240 km/sec), respectively. At u^ = 50 000 ft/sec the correspondence in shock- 
wave angles (a) for a given flow deflection across the shock wave (0) continues to be 
excellent down to very low values of cr and 0 (fig. 32). For u^ = 20 000 ft/ sec 
(6.096 km/sec) this correspondence is not quite so good but is still very satisfactory 
from the normal shock case down to weak shock deflections approaching 20°; for lower 
0's the deviations must increase since at 0 = 0 the shock angle becomes the Mach 
angle which varies in accord with the Mach number of figure 31. 


The density behind the shock wave normalized by the ambient flight density is 
plotted in figure 33. When Uoo = 50 000 ft/sec (15.240 km/sec), the curves for 
T5 = 1000° K and 2000° K follow the curve for T5 = Too within 5 percent from normal 
shock to weak shocks where 0 ~ 20°. As expected the variation is greater when 
Uoo = 20 000 ft/sec (6.096 km/sec), but even here the correspondence is good down to 
the point where oxygen dissociation is essentially complete as indicated by the "bump" on 
the curves at 0-35° for 200 000 feet (60.96 km) and 0 ~ 30° at 300 000 feet 
(91.44 km). After passing across the "bump" the correspondence improves again down 
to 0 ~ 20° where another divergence occurs. 

The post shock temperature (fig. 34) also shows only minor deviations with T5 
down to flow deflections of 20° for u TO = 50 000 ft/ sec (15.240 km/sec). The deviation 
is, as expected, greater for u«> = 20 000 ft/sec (6.096 km/sec) and although it is 
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roughly 500° K for 6 = 20° a wider discrepancy appears for 30° < 6 < 40° where the 
aforementioned dissociation effect is apparent. 

The post shock reciprocal molecular weight ratio is the last shock wave parameter 
illustrated (fig. 35). As in the preceding two figures the correspondence is best at the 
higher velocity and lower altitude. Note the effect of the ambient dissociation Zg - 1 
(fig. 5) which causes the divergence of the weak shock curves as 6 becomes small. 

Reference 19 considered the effects of a constant total enthalpy simulation fairly 
similar to that of the present report. Attention was directed toward larger values of Z 
so that in the approximate energy equation used therein, the thermal energy was negli- 
gible as compared with the dissociation energy (i.e., I15 in eq. (19) was replaced by the 
dissociation energy) and furthermore a hypersonic approximation was invoked requiring 
the kinetic energy normal to a shock wave to be much larger than the ambient thermal 
energy. Nevertheless, the results of reference 19 should still provide some guidelines 
for the present analysis. The symbol oico 0 of reference 19 is equivalent to Z 5 - 1 
herein. Since Zg - 1 is less than 0.04 for equilibrium expansion to region 5, fig- 
ures 13 to 19 of reference 19 show that the effect of predissociation is insignificant for 
the normal- shock condition but is more significant for an attached shock with flow deflec- 
tion of 30°. For this latter condition the hypersonic approximation used in reference 19 
is not on firm ground. Flow -detachment angle and flow deflection at attachment are also 
little changed at Z = 1.04. 

Thus, in regard to simulation under equilibrium conditions one would expect fairly 
good reproduction of the body flow fields with less departure at the higher velocities and 
stronger shock waves associated with blunter bodies. 

Nonequilibrium Case 

Although the energy in vibration or dissociation in the test gas is small compared 
with the flow kinetic energy when the test gas in region 5 is assumed in equilibrium, the 
mere existence of these higher energy levels in the test stream as well as the presence of 
atomic oxygen can appreciably influence the gas chemistry. In reference 20, theoretical 
calculations were made at ambient conditions of T = 294° K, p = 3.15 x 10"8 grams/cm^, 
and Z = 1.0 and 1.1 for shock speeds of 5.9 and 7.2 mm/ psec (u^ = 19 000 and 
24 000 ft/sec). These calculations showed that for the predissociated free- stream 
case the oxygen, nitric oxide, nitrogen, and electron concentrations as well as radia- 
tion from the N2 + and NO* systems were initially increased at the same position 
after the shock when compared with the undissociated free- stream case. However, 
the maximum and equilibrium concentrations were not greatly different, although the 
positions of maximum overshoot and maximum nonequilibrium radiation were closer 
to the shock wave. For a slower shock speed of 4.1 mm/ psec (u,*, = 13 000 ft/ sec) 
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the deviations were more pronounced, especially in regard to the radiation where the 
overshoot was significantly higher when Z = 1.1. 

Although the condition of nonunity in ambient Z was considered in reference 20, 
the case of concurrently varying ambient temperature and velocity was not examined. A 
brief treatment of the latter case may be found in reference 19. No further discussion 
of the finite chemistry and associated nonequilibrium radiation is contained herein since 
generalizations other than the following statements are difficult: (a) The initial ambient 
dissociated oxygen should increase the concentrations closely behind the shock wave of 
oxygen as well as nitrogen and nitric oxide, which require molecular oxygen as an inter- 
mediate product; (b) if the model size is so small that equilibrium is not reached in the 
stagnation region, then the influence of the oxygen atoms in the test section permeate the 
entire flow field; and (c) the deviations will become more significant as the shock wave 
becomes less normal. 


RESUME 

This analysis has considered the effects of simulation wherein blunt-body stagna- 
tion conditions were closely matched for simulated and flight conditions. The simulation 
mode selected was variation of the temperature of the test medium T 5 in the expansion 
tunnel or expansion tube. 

A brief perfect-gas treatment of expansion-tunnel operation with simulation showed 
that both increases and decreases in driver gas pressure P 4 resulted, depending on the 
tunnel area ratio A, driver temperature ratio T^/T^, simulation temperature ratio 
T^/Tx,, and flight Mach number Moo. The use of simulation when was less 

than unity did not appear advantageous from driver-pressure requirements for 
20 < Moo < 50 and 1 < A < 10 ^. 

A more thorough equilibrium real- air study of expansion- tube operation comparing 
simulation at T 5 = 1000° K and 2000° K with exact duplication when T 5 = Too showed 
the following: 

Principal Advantages of Simulation: 

1. The ideal testing time per unit expansion chamber length [^ 2 / 1 S 2 ) increases 
significantly. 

2. The ratio of initial test gas slug length to the expansion- chamber length (zjl 
increases by orders of magnitude. 

3. Equilibrium flow fields and radiation should be similar for the same body shape 
with similarity increasing with increases in velocity and body bluntness. 
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4. Driver-pressure requirements are markedly reduced; or conversely, the per- 
formance potential in terms of flight velocity and altitude for a given driver are greatly 
increased. 

5. The efficiency of an arc-heated helium driver is increased. 

6. The minimum pressure required in the expansion chamber is increased an order 
of magnitude and thus vacuum pumping requirements are reduced. 

Principal Disadvantages of Simulation; 

1. The reciprocal test- section molecular weight ratio Z5 departs from unity for 
Tg = 2000° K. The resulting degree of dissociation (Z5 - l) can contribute up to 3 per- 
cent of the total enthalpy (H5) for the worst case evaluated, namely, flight velocity 

Uoo = 20 000 ft/sec (6.1 km/sec), altitude = 300 000 ft (91.44 km). 

2. The higher initial thermal and (for Z5 ^ l) chemical energy of the test- section 
gas for T 5 = 1000° K and 2000° K may significantly influence finite rate chemistry. 
Consequently, nonequilibrium flow fields and radiation may vary with simulation; the 
variation will probably be greater for lower velocities, lower shock wave angles, lower 
densities, and smaller models. 

3. The degree of dissociation behind the primary shock wave (Z2 - l) increases 
significantly with simulation. If freezing does not occur in the unsteady expansion, then 
large Z2 is not important. However, if appreciable freezing occurs, then Z5 would 
be increased above the equilibrium values. 

4. Larger lengths of intermediate chamber (l and driver chamber are 
required for a given l gg. However, for high values of Uoo where the advantages of 
simulation become marked, the lengths Z j) and Z sj are still approximately two orders 
of magnitude less than Z 32* 

5. The test- section Mach number is reduced, but this reduction is not a significant 
penalty as long as flow is still hypersonic. 

6. The test- section viscosity is increased. The detrimental extent of this effect on 
test core size and attenuation is not known. 

The relative significance of the aforementioned advantages and disadvantages will 
vary with the facility size and equivalent free- stream ambient conditions desired. Con- 
sequently, no clear-cut general conclusion regarding the use of simulation is possible, 
although the advantages of simulation increase as free-stream velocity or density (or 
both) increase. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., January 21, 1966. 
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Figure 5.- Reciprocal molecular wight ratio Z 5 as ^function of altitude. 





























Altitude, ft (km) 


— 100 000 
- 150 000 
—200 000 
-250 000 
-300 000 


(30. 48) 
(45. 72) 
(60.96) 
(76.20) 
(91. 44) | 


40 50 



Flight velocity, u^, 10 


(b) T 5 = 1000° K 
Figure 11.- Variation of p 2 wi 




20 30 40 51 


(c) T 5 = 2000° K. 


flight velocity. 






Reciprocal molecular weight ratio in region 2, 



Figure 12.- Variation of reciprocal molecular weight ratio in region 2 with flight velocity. 
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Figure 13.- Variation of shock velocity with flight velocity. 
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Figure 30.- Effect of simulation on maximum velocity capability of expansion tube using three typical driver modes. 
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Figure 34.- Variation of static temperature behind shock wave with flow deflection angle across shock wave. 
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Figure 35.- Variation of reciprocal molecular weight ratio behind shock wave with flow deflection angle across shock wave. 








